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I. INTRODUCTION 

Systems with coupled scalar fields hold their certain interest in physical applications. In particular, such systems 
are used in constructing particle models and their interactions within the framework of quantum field theory In 
different aspects such systems have been already repeatedly considered, see for example Rcfs. Here we present 

a qualitative study of one of such systems having the potential energy in the form: 

= ^(^2 - + ^(^^ - + + const, (1) 

where ip,ip are two real scalar fields (usual or phantom/ghost ones), Ai, A2, A3, mi, m2 are some constants. This 
potential had been used earlier by us in treatments of models of cosmological and astrophysical objects in general 
relativity. These researches showed that: (a) for the four-dimensional case there exist regular spherically and cylin- 
drically symmetric solutions jol-fTH. and also cosmological solutions [13, [13] both for usual and phantom scalar fields; 
(b) for the higher dimensional cases there exist the thick brane solutions fl^-flTj supported by usual and phantom 
scalar fields. 

From the physical point of view, these solutions exist because of the special form of the interaction potential ([T]) 
having two local and two global minima. It means that there are two different vacua. At the infinity, as the radial 
coordinate x — >■ ±00, these scalar fields are located in that vacuum in which their are in the local minimum. The 
existence of regular solutions with finite energy is only possible for certain (eigen) values of parameters of a problem. 
Such eig envalue problems were solved by us using the shooting method (see for details of the shooting procedure in 
Ref. [nl). 

Note that the type ([T|) potential has been also used in the paper [Tsj for modeling superconductivity using two 
coupled Ginzburg-Landau equations. When the interaction between the fields in ([1]) is excluded, i.e. at A3 = 0, one 
has two uncoupled type Ginzburg-Landau equations. On the other hand with an account taken of the interaction 
there are new effects, in particular, the presence of regular solutions, not present in the case of one Ginzburg-Landau 
equation. 

Use of spherical and cylindrical symmetries in the above papers leads to obtaining the sets of coupled non- 
autonomous ordinary differential equations whose qualitative study is quite complicated. Here we consider a simplified 
problem when the system with the potential ([T|) is examined in cartesian coordinates. This allows to get an autonomous 
system of two second order ordinary differential equations for which it is possible to perform the qualitative analysis 
and estimate the general behavior of the system. 
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II. QUALITATIVE ANALYSIS 

Let us consider the physical system with the potential ^ whose Lagrangian can be presented in the form 



where ei,e2 = il, and the plus sign refers to usual scalar fields, and the minus sign - to phantom/ghost ones. The 
energy- momentum tensor of the system is: 

= eid.ipd'^ip + esa.^/'^V - 5'IL. (3) 
The corresponding field equations in cartesian coordinates can be written as: 

= , (4) 

ei dip 

Choosing the initial value of the scalar field Lp{r — Q) — ip^ and introducing the dimensionless variables (j) = ip/ipo, x — 

ip/cpojX = \/A^Lpor, = mi/cpQ, /i2 = 1122/^0, and also Ai — Ai/A3,A2 = A2/A3, let us rewrite the system in the 
form (taking into account the expression for the potential from ([T])): 

0' = ^, (6) 

X' = V, (7) 

z'^-cb[x' + Xi{cb'~f^l)], (8) 

-X[0^ + A2(x^-M^)]. (9) 
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The fixed points of this system are: 





{2; ^ 0, in- 0, X 


^0,0 


Ml}, 


B ^ 






-Ml}, 


C = 


{z 0,v^ {),4> 




M2}, 


D = 


{2 -> 0, u -> 0, (/) 


->o,x 


-> -M2} , 


E = 


{z -> 0,w ^ 0,0 


^0,x 




F = 


|z ^ 0,u ^ 0,<^ 




'A2M2 - A1A2M? 


1-A1A2 



(10) 

(11) 

(12) 

(13) 
(14) 



,x-W^4^r (15) 
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The points A, B refer to local minima, the points C, D are the global minima, E is the local maximum, and the points 
F refer to saddle points. Let us designate the values of the potential ([1]) at these points as Vi where the index i 
corresponds to letters from A to F, and Va = Vb and Vc = Vd- 
Obviously that 

max{VA:Vc)<VE. (16) 

It is also possible to find the following relations: 

Ai(-M? + A2/i2)^ 



Vf-Va^ 
Vf-Vc^ 
Vf-Ve = 



4(l-AiA2) ' 

A2(-Mi + Ai/ii)^ 
4(l-AiA2) ' 

A1A2 [mi(Aimi - ^4) + ^4{X2^4 - mi)] 

4(l-AiA2) 
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Dividing the third expression by the first one and the second one and taking into account the condition (jl6p . one can 
show that the conditions guaranteeing the existence of the local and global minima are: (i) for the local minimum; 
Ai > 0,^1 > A2/I2; (ii) for the global minimum: A2 > 0,^| > Xifi'f. 

Next, to determine the type of the fixed points it is necessary to make a characteristic matrix of the system ([S])-®- 
Using this matrix, one can write out the characteristic forth order algebraic equation. Then the values of roots of this 
equation kj determine the type of the fixed points. In our case we have the following roots: 

(i) at the points A, B (the local minimum): 

2 1 
fci = —Xinl, ^2 = — (/^i - A2/^2) : ^3,4 = 1- 

(ii) at the points C\D (the global minimum): 

1 2 

fci = — (^2 - ^11^4) 7 ^2 = — A2/i2, ^3,4 = 1. 

ei £2 

(iii) at the point E (the local maximum): 

fcl = - — Xiflj, k2 = - — A2/i2, ksA = 1- 

ei £2 

(iv) at the points F (the saddle points): 

fcl. 2 = 7t\ -t| - e2Mi^i^2 - eifJ-jXiXl + eififXiX-z + £2/^2^1 A2 

eie2(AiA2 - 1) 1^ 

±<^XiX2 {4eie2(AiiA2 - - Xi^il){XiX2 - 1) + A1A2 K(£2Ai - ei) + ^^^£1^2 - £2)]^} }, 

^3,4 = 1. 

Taking into account the above mentioned conditions of existence of the local and global minima, one can make the 
classification of the fixed points presented in table ID 

TABLE I: The classification of the fixed points of the system ([B])-® for different values of the parameters ei and £2- 





Points A, B 


Points C, D 


Point E 


Point F 


ei = 1, £2 = 1 


Unstable node 


Unstable node 


Saddle 


Depends on the values of 
the parameters Ai, A2 


ei = 1, £2 = —1 


Saddle 


Saddle 


Saddle 


Depends on the values of 
the parameters Ai, A2 


£1 — —1, £2 = 1 


Saddle 


Saddle 


Saddle 


Depends on the values of 
the parameters Ai, A2 


£1 = —1, £2 = —1 


Saddle 


Saddle 


Unstable node 


Depends on the values of 
the parameters Ai, A2 



III. NUMERICAL ANALYSIS 

From the point of view of obtaining a set of solutions (but not only of two integral curves as in a case of saddle fixed 
points) type "unstable node" fixed points are more interesting. In the model under consideration, they are situated: 
(i) at the points of the local {A,B) and global {C,D) minima at positive £i,£2, i-c. for usual scalar fields; (ii) at the 
point of the local maximum E at negative £1, £2, i-e. for the case of phantom/ghost fields. Note that in the latter case 
negative ei, £2 effectively correspond to the system with usual fields but with a reversed sign of the potential ([1]). In 
this case the point of the local maximum E becomes the point of the local minimum, and solutions asymptotically 
tend to that point. 

All the fixed points, being the stationary points of the system (jH])-®, arc situated at a; = ±00. Then static 
solutions, if they exist, should start from these points. One of physically interesting problems is a Z2 symmetric 
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TABLE II: The initial values of xo and the corresponding values of the parameters ^i,/X2 for the system ((Gjl-®. 



# 


4>o 


xo 


Ml 


M2 


M 


1 


1.0 


0.3 


1.25104535 


1.1056305 


0.0441854 


2 


i.U 


V U.z 


i.4o44oo 1 


1 1 o^onf;o 
i.io/o9DO 


U. 148/09 


3 


1.0 




1.736266 


1.30665 


0.41906 


4 


1.0 


\/a6 


1.9628773 


1.40650056 


0.76536 


5 


1.0 


^/a8 


2.158048 


1.495301394 


1.17074 


6 


1.0 


1.0 


2.33213652 


1.5764432135 


1.62578 


7 


1.0 


^L2 


2.4908109 


1.6518053896 


2.1246 


8 


1.0 


\/L4 


2.63757479 


1.7225756427 


2.66294 




FIG. 1: The scalar fields 0, x from the system (|6])-(|9]) for the FIG. 2: The dimensionless energy density e of the system (|6}- 

different initial values of xo taken from table [III Asymptoti- ([9]) from (|17[l for the different initial values of xo taken from 

cally, as a; — >■ ±cxj, the scalar field x — >■ 0, and 4> goes to values table [Til The top line corresponds to the greatest xo ~ Vl.4, 

of /ii from table|II]corresponding to the local minimum A from and the bottom one - to the lowest xo = 0.3. 



solution. In this case the symmetry plane is chosen at x = where the derivatives ^'(0) = x'(0) = 0. Examples of 
such solutions are presented in Fig. [1] Using the energy-momentum tensor ([3]) and the above dimensionless variables, 
the dimensionless energy density can be derived in the following form 




where the constant const from the potential ([T]) is chosen equal to — (A2/4)/i| to make the energy density to be equal 
to zero at infinity. Using this expression, one can plot the corresponding graphs for the energy presented in Fig. [2] 
The total energy (mass) of the system is defined by the expression 

/oo 
e{x)dx. 
'OO 

Calculating this integral, the values of the total energy presented in the last column of table HI] have been found. Also, 
using the table, one can plot the graphs of dependence of the parameters /ii, fi2 and the total energy M on the initial 
values of Xo presented in Fig. |31 The corresponding phase portraits are shown in Figs. I4I5I 

The system (H])-© allows introducing another dimensionless variables. Namely, introducing the dimensionless 
variables cj) = ip/mi,x = ip/mi,x = ^/A^mir, fi = m2/mi, and also Ai — Ai/A3,A2 = A2/A3, one can rewrite this 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



FIG. 3: The dependence of the parameters of the system /ii,/i2 and the total energy M on the initial values of xo for the 
system ©-([U}. The data are taken from table HIl 




FIG. 4: The phase portrait of the system ©-([U for the scalar 
field (j) plotted using the parameters from table [IT] The points 
Ai correspond to the fixed point A from (jlOp . All solutions go 
to these points at a; — >■ ±cxj starting with the different initial 
values at the point O corresponding to the origin of coordinates 
1 = 0. The index i runs over i = 1..8 in accordance with 
numeration from table HIl 




FIG. 5: The phase portrait of the system (|6}-([9]) for the scalar 
field X plotted using the parameters from table [IT] The points 
Oi correspond to the origin of coordinates x = Q with the 
different initial values taken from table [Til Asymptotically, as 
X — >■ ±cxD, the solutions tend to the fixed point A from (|10p 
where x(±oo) = 0. The index i runs over i — 1..8 in accordance 
with numeration from table HIl 



system as follows: 



0' 


= z, 


x' 


= V, 


-1 


1 


z 






ei 




1 


v' 






£2 



4> [t + Ai _ 1)] , 



(18) 
(19) 

(20) 
(21) 



The values of the parameter /i and initial conditions 0(0) and x(0) at which regular solutions do exist can be obtained 
from the values presented in table HIl by corresponding rescaling the variables: 0o = ^o/mi , Xo = Xo/a*i , M = ^J■2/^J■l,x = 
/UiX, M = M/^\. The corresponding new values of the mentioned parameters are presented in table HIB Using this 
table, one can plot the graphs of dependence of the initial values 0oi Xo on M presented in Fig. |5| 

Summarizing, here we have considered the system with two non-gravitating coupled scalar fields in cartesian coor- 
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TABLE III: The initial values 0o, Xo and corresponding to them values of the parameter ^ for the system (|18p - (|21|l . 





XO 




M 


0.799332 


0.239799 


0.883765 


0.0225663 


0.687528 


0.307472 


0.816713 


0.0481829 


0.575949 


0.364262 


0.752563 


0.0800622 


0.509456 


0.394623 


0.71655 


0.101201 


0.463382 


0.414461 


0.692895 


0.116487 


0.428791 


0.428791 


0.675965 


0.128174 


0.401476 


0.439795 


0.66316 


0.137485 


0.379136 


0.4486 


0.653091 


0.145127 



0.8 



0.6 



0.4 



0.2 



0.0 

0.65 0.70 0.75 0.80 0.85 V- 

FIG. 6: The dependence of the initial values cj>o,xo and the total energy M on the value of the parameter p for the system 
dlSJ-lETJ- The data are taken from table Hill 

dinates. For such a system, as well as in general relativity, the task of finding regular solutions amounts to searching 
eigenvalues of the parameters of the model. The model contains six available parameters: two initial values of the 
scalar fields (po, V'o and four free parameters Ai, A2, mi, 7712 (A3 can be always excluded by redefinition of the parame- 
ters Ai, A2, and the initial values of derivatives <p'(0), V''(0) are chosen to be equal to zero for obtaining Z2 symmetric 
solutions). Then for obtaining regular solutions it is necessary to find eigenvalues of only two of these six parameters. 
For example, we have been sought the eigenvalues of the parameters jii and 112 m the system (jH])-®. For the obtained 
eigenvalues regular solutions start at a; = and tend to the fixed point A corresponding to the local minimum of 
the system (for the case of usual scalar fields considered here, see Figs. [1] U and [5]) , and to the local maximum (for 
phantom/ghost fields). One can see from Fig. [3]that there exist some lowest eigenvalues of the parameters /ii and ^2 
at which the total energy of the system M goes to zero. This corresponds to the existence of some critical /ii and ^2 
at which physically sensible solutions with a nonzero total energy still exist. Similarly, for the system ([T5| - (PT|) there 
exists some critical eigenvalue of the parameter ^ at which M — >■ as well (see Fig. |6|). 

Type "unstable node" fixed points allow the existence of sets of solutions starting from these points at x = ±00 
both for usual fields (the fixed points A,B and C,D) and for phantom/ghost scalar fields (the point E). Thus the 
qualitative analysis shows that from the point of view of a possibility of obtaining regular localized solutions the 
system with two coupled scalar fields in question seems to be quite perspective. The previous studies from Refs. @- 
[17| show that inclusion of gravitational fields does not change the qualitative behavior of solutions. Then one would 
expect, for example, that in the presence of gravitational fields the lower restriction on the values of the parameters 
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/ii and /i2 giving physically sensible solutions will also exist. 
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